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Abstract 



The surface free energy of ideal hard rods near curved hard surfaces is 
determined to second order in curvature for surfaces of general shape. In 
accordance with previous results for spherical and cylindrical surfaces it is 
found that this quantity is non-analytical when one of the principal curvatures 
changes signs. This prohibits writing it in the common Helfrich form. It is 
shown that the non-analytical terms are the same for any aspect ratio of 
the rods. These results are used to find the equilibrium shape of vesicles 
immersed in solutions of rod-like (colloidal) particles. The presence of the 
particles induces a change in the equilibrium shape and to a shift of the 
prolate-oblate transition in the vesicle phase diagram, which are calculated 
within the framework of the spontaneous curvature model. As a consequence 
of the special form of the energy contribution due to the rods these changes 
cannot be accounted for by a simple rescaling of the elastic constants of the 
vesicle as for solutions of spherical colloids or polymers. 



PACS numbers: 82.70.-y, 68.10. Cr, 87.22.Bt, 62.20.Dc 
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I. INTRODUCTION 



Recently Yaman et al. |I]J3] obtained the surprising result that the surface tension 7 of 
a fluid of thin hard rods at a spherical or cylindrical wall of radius R is non-analytical at 
curvature c = 1/R — 0. They found different expressions for 7(c) for positive and nega- 
tive curvatures. This analysis that assumed ideal particles has been extended by Groh and 
Dietrich || by taking into account the interactions between the particles. This leads to a 
substantial change of the curvature dependence at medium and high particle densities, but 
a small singular contribution probably still remains. In the present work we return to the 
ideal limit but generalize the problem in two respects. First, we explicitly study arbitrarily 
shaped surfaces with locally varying curvatures. This is especially important because due to 
the mentioned singularity the curvature dependence of the surface free energy cannot have 
the common Helfrich form |4j]. Therefore it is not possible to derive the general expression 
from the special cases of spherical and cylindrical surfaces. Indeed, we find different analyt- 
ical expressions depending on the signs of the local curvatures which take on an unexpected 
form especially at hyperbolic points. Second, we allow for a finite thickness of the rod-like 
particles, with the known result for spherical particles included as a limiting case. 

One might ask if and how the curvature dependent surface free energy is experimentally 
accessible. This is indeed the case for surfaces that are free to adjust their curvature to the 
conditions provided by the surrounding liquid. Such surfaces are given by vesicles which 
consist of closed liquid bilayer membranes in an aqueous solution. Their shapes are experi- 
mentally observable by phase-contrast microscopy |J. Theoretically vesicles are commonly 
modeled as two-dimensional continuum surfaces with bending elasticity. In the simplest 
model the elastic energy of the membrane has just the Helfrich form. Different classes of 
equilibrium shapes have been determined as local energy minima (for some especially in- 
teresting examples see Ref. ||). For the most relevant parameter ranges the corresponding 
phase diagram has been worked out in detail J7| . 

The presence of other particles in the solution gives rise to an additional energy con- 
tribution. Due to a change of the effective spontaneous curvature the vesicle may tend to 
bend towards or away from the solute particles. The effects of absorbed and free polymers 
and of spherical colloidal particles have been reviewed by Lipowsky et al. J8|. Here we study 
the corresponding problem for rod-like particles. In contrast to the other cases, here the 
total elastic energy including the solute effects no longer has the Helfrich form. Thus new 
equilibrium shapes arise and the phase diagram is modified in a way that cannot be reduced 
to a simple renormalization of the parameters. 

II. SURFACE FREE ENERGY AT ARBITRARILY SHAPED WALLS 

We study a system of hard sphero cylinders of length L and diameter D in the presence of 
a hard wall of general shape whose principal radii of curvature R\ and R2 are large compared 
to L and D. The interactions between the rod-like particles are neglected, i.e., we consider 
the dilute limit. In this section we calculate the surface free energy of the rods, which is then 
used in the second part of the paper to determine equilibrium shapes of vesicles immersed 
in the rod solution. 
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Using the grand-canonical density-functional it can easily be shown 0] that the surface 
contribution Q s to the grand-canonical potential of the fluid is given by 

where p(r, uj) denotes the density of rods whose center of mass is at the point r and whose 
orientation is uj = (8, <fi) and pb is the particle density of the bulk fluid far away from the 
wall. Due to the hard wall potential n(r, uj) : = 4np(r, uj)/p b has either the value or 1, 
depending on whether the configuration (r, uj) is forbidden by the presence of the surface or 
not. The ideal particles we assumed only "feel" the wall within a layer of thickness (L + D)/2 
so that the spatial integration in Eq. (|I|) can be restricted to this layer. We introduce a local 
coordinate system at each point of the surface whose z axis points along the normal direction 
into the rod solution and whose x and y axes are aligned with the principal directions. As 
we are only interested in the leading terms in the curvature we approximate the position of 
the surface by 

*~KS + £) 

where we introduced the convention that a curvature q = 1/Ri is positive if the surface 
curves away from the fluid. Now Q s can be written as a surface integral over the free energy 
density uj s : 

Q s = J dSuj s (c 1 ,c 2 ) (3) 

with 

ML+D)/2 1 , 

(3uj s {c 1 ,c 2 ) = pb I dzJ(z,c 1 ,c 2 )— duj{l-n{z,ci,c 2 ,uj)). (4) 







The Jacobian 



t/ \ (z + Ri)(z + R 2 ) , , 2 

J{z, ci, c 2 ) = i g± ^ = 1 + (ci + c 2 )z + Cl c 2 z 2 (5) 

ti\ti 2 

takes into account the change of the tangential area element dS with the normal distance 
z. The range of orientations uj for which the particle intersects the surface (i.e, for which 
n(z,ci,c 2 ,uj) = 0) for given curvatures c\ and c 2 and normal distance z can be found by 
straightforward geometry. 



A. Infinitely thin rods 

We first consider the limit of infinitely thin rods (D = 0). At given azimuthal angle 
around the normal direction we determine the maximum allowed value x max for the cosine of 
the polar angle, at which the rod just touches the surface. This is a two-dimensional problem 
within the plane = const that intersects the surface in the parabola z s = —cp 2 (p is the 
tangential coordinate within the plane) with 
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c = c(4>) = c\ sin <p + c 2 cos <p (6) 

Note that c(0) can be positive or negative depending on the signs of c\ and c 2 . Figure [I] 
shows that for positive curvatures there are two possibilities how the thin rod can touch the 
surface: tangentially for z < z c and with its end for z > z c . It is easy to show that 

^W/ 2 = ^/2^+0(c 3 / 2 ), 



^ " { ~& + S + S = T + (f - (£ - f )c 2 + 0( C 3), * > z c (7) 



with 



cL 
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— (-1 + Vl + c 2 L 2 ) = — + 0(c 3 ) (8) 
4c 8 



For negative curvatures the rod cannot approach the surface closer than z c = —cL 2 /8 so 
that x max = for z < z c while for z > z c one has 

* max = 2 { + - j)c + (j ~ f )c 2 + 0(c 3 ) (9) 

as before. Note that x max is analytic at c = only if z > z c while different expressions apply 
for c 5s if z < z c . The surface free energy density now follows by integration over 0: 

— -(ci,c 2 ) = — / d(j) dzJ(z,c u c 2 )(l - x max (z,c((f)))). (10) 
Pb ^ Jo Jo 

Up to terms quadratic in the curvature one obtains 

(3uj s L 



Pb 4 

and 



for ci, c 2 > (11) 



(3u s _ L L 6 n _ t x2 4 
~ 4 ~ 128 



((ci + c 2 ) - -cic 2 ) for ci, c 2 < (12) 



If Ci > and c 2 < the function c(0) changes sign at the angle O — arctan a/— c 2 /ci and 
the integration must be split accordingly. After some algebra one finds 



f3uo s L L 
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3 v / -cic 2 (ci + c 2 ) + (3c 2 + 2cic 2 + 3c 2 ) arctan \J-c 2 /ci (13) 



p 6 4 192vr 

for ci > 0, c 2 < 0. 

Of course the same formula with c\ and c 2 interchanged applies for c\ < and c 2 > 0. 
Thus depending on the signs of the principal curvatures u s is given by Eq. (|ll|), (0), or 
(|j~3|). Equation ([3~2| ) was already derived by Yaman et al. (see Eq. (A. 9) in Ref. 0) who 
also showed that Eq. (|TTD is exact for a convex surface to all orders in the curvature. Our 
new result Eq. ( |T3"D provides a continuous connection between Eq. ([TTJ) and Eq. (TT2]). This is 
demonstrated in Fig. |2| where u s is plotted for c\ = ccosip and c 2 = c sin ip as a function of 
ip for fixed c. Obviously the presence of the rods favors surfaces with two negative principal 
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curvatures. Terms linear in the curvature are absent in all cases. The quadratic terms are 
non-analytic if one of the q's changes sign, which, inter alia, precludes writing the curvature 
dependence of the surface free energy in the common Helfrich form 

Fuel = J dS{^K( Cl + c 2 - C ) 2 + R Cl c 2 ) (14) 

with the bending rigidities k and R and the spontaneous curvature Cq. If one restricts 
oneself to equal signs of C\ and c 2 we have rigorously shown here that the Helfrich form can 
be applied with /3n/pt = — L 3 /64 and (3R/pb = L 3 /96 for negative and, trivially, k = R = 
for positive curvatures. However, surfaces of vesicles typically exhibit regions with differing 
signs, so that a global mapping to the Helfrich form is no longer possible. 



B. Rods of finite thickness 

The case of spherocylindrical particles with finite thickness D can be reduced to the 
problem for D = discussed above if one realizes that the surface of the spherocylinder 
consists of all points that have a distance D/2 from the line segment of length L which 
connects the centers of the hemispherical caps. Thus if the spherocylinder touches the surface 
S this line segment touches a parallel surface S' shifted by D/2. The curvature radii of 5" 
and S are simply related by R[ = Ri + D/2, i.e., c[ = q/(1 + CiD/2) |J, so that for z > D/2 

n(z, ci, c 2 , D) = n(z - D/2, c x /{\ + c x D/2), c 2 /(l + c 2 D/2), 0). (15) 

The particles cannot approach the wall closer than z = D/2. The inaccessible range yields 
the L independent contribution 



(1) pD/2 
s 



f u/z D D 2 D 6 

/ dz J(z,d,c 2 ) = — + — (ci + c 2 ) +— cic 2 (16) 

Pb Jo Z 8 24 

to the free energy density. By using the fact that 

J(z' + D/2, ci, c 2 ) = J{z', c[, c 2 )(l + Cl D/2){l + c 2 D/2) (17) 
one obtains for the remaining z integral from D/2 to (L + D)/2 

u>f\ Cl ,c 2 ,D)=(l + C ^ (1 + ^)^,40). (18) 

Now the above results for infinitely thin rods can be employed to give the final result 
(3lo s D L D 2 DL., , ,D 3 D 2 L 

— = — + -+ — + (ci + c 2 + — + cic 2 + A/ a 19 

p b 2 4 v 8 8 A ; V 24 16 ' J K ' 

with 

Af s = 0, d,C2>0 (20) 
L 3 4 

Af s = - — (( Cl + c 2 ) 2 -- Cl c 2 ), Cl ,c 2 <0 (21) 



L 
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1927T 



3 V - c ic 2 (ci + c 2 ) + (3c 2 + 2ciC 2 + 3c 2 ) arctan — - 



ci > 0,c 2 < 0. (22) 
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It is remarkable that all the additional terms which contain powers of D are analytic in 
c\ and C2- The singular contributions A/ s are exactly the same as in the limit D — > 0. In 
this sense no qualitatively new features arise when the thickness of the rods is taken into 
account. Therefore we only consider thin rods in the second part of this paper. The special 
case of hard non- interacting spheres is also contained in the above analysis for L = 0. Only 
in this case the full uj s is analytic and thus lends itself to the Helfrich expansion. Our result 
for the spheres is in agreement with previous calculations JT|,[TTJ | . 



III. VESICLE SHAPES IN A ROD SOLUTION 

We now consider a vesicle immmersed in a solution of rodlike colloids. It is assumed that 
the rods are not absorbed on the membrane and that their interaction can be approximated 
by the hard wall model. The bending energy of the vesicle consists of two contributions: 
First, the internal bending rigidity of the membrane which is described in the so-called 
spontaneous curvature model as given by Eq. (|i4|) (For other possibilities see, e.g., |5|,p!T|]). 
Second, the presence of the rodlike particles gives rise to the surface free energy determined 
in the previous section. If the vesicle is immersed in a solution of polymers [12| or diluted 



spherical colloids this second contribution will have the same analytical form as the first 
and therefore lead to a renormalization of the rigidity coefficients and the spontaneous 
curvature of the "free" vesicle. Then the equilibrium shapes are those determined for the 
pure spontaneous curvature model |7j and the addition of the solute only shifts the system to 
a different point in the phase diagram. However, in the present case the energy contribution 
due to the rods cannot be written in the Helfrich form so that new equilibrium shapes will 
arise. As shown in Sec. |IIB| a non-zero thickness of the rods can also be taken into account by 
using renormalized coefficients in Fuel- Therefore in the following we only use the expressions 
for infinitely thin rods. 

Because the energy scale associated with the bending rigidity k is much smaller than the 
energies necessary to change significantly the area A or the volume V of the vesicles, these 
quantities are regarded as fixed. Thus the equilibrium shape is determined by minimization 
of F tot = F He i + Q s under the constraint of fixed A and V. Due to the scale invariance of 
F to t the equilibrium shape depends only on the dimensionless parameters 



V „ A Pb L 



v = ^a^> c ° = c °y^ x= ^k~- (23) 

A sphere is characterized by v = 1, while for any other shape v < 1. The parameter x 
measures the relative importance of the contribution Q s due to the rods. The Gaussian 
curvature term proportional to R in Eq. (|TJ]) is constant for topologically equivalent shapes 
and will be omitted since we consider only shapes with the same topology as the sphere. We 
will focus on the region v < 1 and moderate values of Cq where the equilibrium shapes of 
the pure spontaneous curvature model are axisymmetric prolates or oblates with up-down 
symmetry, separated by a first-order transition 0. 

Rather than deriving and solving the exact Euler-Lagrange shape equations we 
parametrize the surfaces by simple two-parameter functions. In cylindrical coordinates 
(r, z, 0) around the axis of symmetry we write |TJ[ 
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r(z) 
z(r) 



a\/l — z 2 \J\ + bz 2 prolates 
ayl — r 2 Vl + 6r 2 oblates 



(24) 
(25) 



These shapes reduce to ellipsoids for 6 = while increasing the value of b leads to a bulging 
which finally yields non-convex shapes for b > 1. The parameter a is the aspect ratio of 
the shape, i.e., the ratio of its diameter in the equatorial plane and its length along the 
symmetry axis (or vice versa for oblates). We emphasize that the bending energy does not 
depend on the size of the vesicles so that the length of one axis can be arbitrarily set to 
unity. The corresponding vesicle volumes are 



V, 



pro 



Vr 



obi 



Air 

T 

2b 



a 2 (l + 



(b + l) 
2Vb 



vr 6-1 

2 +arcsin 6TT 



and the areas 



.4 



pro 



A 



obi 



Ana 



An 



dzy/(l - z 2 )(l + bz 2 ) + a¥(l - 6 + 26^ 2 ) 2 



drr\ 1 + a 2 r 2 



(1 -6 + 26r 2 ) 2 
(1 -r 2 )(l + 6r 2 ) 



(26) 
(27) 

(28) 
(29) 



The principal curvatures are 

<$ ro = r(z)y/l + r'(z)*, 



pro 
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for prolates and 



Mbl 



z'(r) 



„obl 



{l + r'(z) 2 ) 3 / 2 
? i (z) 

(1 + z\rffl 2 
z"{r) 



(30) 



(31) 



for oblates. All curvatures are positive for 6 < 1. In the prolate case c^ ro becomes negative 
for z < z c where z c is the positive solution of 

- 1 + 6V(-3 + 2z 2 ) + 6(1 - Qz 2 + 3z 4 ) = 0. (32) 

For oblates three ranges must be distinguished: for r < r cl , where r cl is again determined by 
Eq. (|3"2"D, both curvatures are negative. Between r cl and r c2 = a/(6 — l)/(26) c ^ 1 is positive 
and c° w is negative, while for r > r c2 both curvatures are positive. The rod contribution to 
the elastic energy is computed from 



dS io s (ci, c 2 ) 



Att J dzr(z)y^l + r'(z) 2 u> s ((^ ro (z), (^°{z)) prolates 
4tt /q 1 dr r y/l + z'(r) 2 uo s (cf (r ) , cf (r ) ) oblates 



(33) 



The equilibrium shapes follow by numerical minimization of the total elastic energy 
Ftot{a{b,v),b) with respect to 6 where a(b,v) is obtained by numerical solution of the first 
equation of Eq. (p3|). 
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If the rods are inside instead of outside the vesicle the signs of the curvatures q have to 
be reversed in the computation of Q s . For the case of particles on both sides of the membrane 
^s(ci,c 2 ) must be replaced by cj s (ci,c 2 ) + uj s (— ci, — c 2 ) which, interestingly, has the form 
given by Eq. ([12]) for all signs of the q. Hence here the Helfrich form is valid, as has been 
surmised by Yaman et al. so that for Co = the equilibrium shape does not depend on 
the rod concentration x. 

We first discuss the results for zero spontaneous curvature. If b > 1 the shape can be 
characterized by a "bulging" parameter t — (b + 1) /26 1//2 which is equal to the ratio of the 
maximum of r(z) [z(r)] and its value in the equatorial plane [along the symmetry axis] for 
prolates [oblates] . In Fig. [5] the quantities a and t are shown as a function of x for a prolate 
and an oblate at v values in the vicinity of the phase transition. Rods outside a prolate tend 
to decrease a and increase t thereby narrowing the waist of the vesicle. Similarly an oblate 
develops stronger "dips" at the symmetry axis in order to increase the range and degree of 
negative curvature. These results confirm our previous observation that the surface prefers 
to bend towards the rods. Figure f| shows some examples for shapes without and with rods 
outside and inside the vesicle. The changes of a and t are much smaller for prolates (0.4-1.3% 
at x — 30) than for oblates (10-15%). 

Identifying the most probable shape with the lowest energy shape is problematic when 
the energy is a non-analytic function of the shape which might be non-quadratic for small 



deviations from the energy minimum [ Hf| . However, numerically the function Q s (a(b, v), b) 



shows no traces of non-analyticity. In fact, the integrand oj s (ci, c 2 ) is singular only on a set 
of measure zero (on the curves determined by cic 2 = 0) so that the singularity is probably 
removed by the integration. 

Without rods the phase transition between prolates and oblates takes place at v = 0.648 
within our approximation which is very close to the value v = 0.651 obtained by an exact 
minimization |7j . The coexisting shapes (see Fig. |j) both are strongly non-spherical with an 
aspect ratio a ~ 0.17 which indicates a pronounced first-order character of the transition. 
As shown in Fig. |5| the rods shift the transition to larger v if they are outside the vesicle and 
to smaller v inside. As explained above rods on both sides effectively just change the value 
of k but do not influence the phase diagram. 

We now turn to the more general case of non-zero spontaneous curvature Cq. Negative 
values of Cq favor oblate shapes so that the transition point moves to higher volume to 
surface ratios v (Fig. ||). It has been shown rigorously that the phase boundary approaches 
v = 1 at c = — 6/5 [0,0- By a series expansion of the elastic energy around v — 1 one 
finds that this also hold within the present parametrization. The presence of rods outside 
the vesicle shifts the transition line to larger v. However, no shift occurs for v > 0.87 because 
in this range both coexisting shapes are convex (b < 1) so that Q s = const. Rods on both 
side of the membrane give rise to an additional Helfrich like term with (3K ro d = — pbL 3 /64 
and c 0t r O d = 0. Therefore the effective bending rigidity is k' = k + K rod and the effective 
spontaneous curvature is 

(34) 



l-x/64 

For this reason the curve for particles on both sides differs from the curve for x = in 
Fig. |6| (except at Co = 0) and reaches v = 1 at cq = —6/5(1 — x/64), whereas both curves 
are identical if plotted as a function of c' instead of cq (Fig. 0). We note that k' becomes 
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negative for x > 64. In this case the present analysis breaks down and higher order terms 
in the curvature must be taken into account. If the rods are restricted to the inside of the 
vesicle the same effective Helfrich form applies as long as only convex shapes are considered. 
This explains why the corresponding phase boundary is equal to that for rods on both sides 
at large v. However, for smaller v shapes with differing signs of the curvatures occur which 
cannot be accounted for by a simple rescaling of the Helfrich coefficients and which shift the 
transition line in Fig. |7| to lower v compared to the pure Helfrich case. 



IV. DISCUSSION 

In summary, we have shown that if rod-like particles are present on the outer or inner 
side of a vesicle its equilibrium shape changes. Because the curvature dependence of the 
surface free energy of the solutes cannot be written in the Helfrich form their effect cannot 
be described by a simple rescaling of the bending rigidity coefficients. We have quantitatively 
computed the shape changes and the shift of the prolate-oblate transition line in the phase 
diagram. 

In order to decide whether these effects are large enough to be observable in experiments 
an estimate for the quantity x = p^L? / (3k is needed. The bending rigidity is typically of 
the order of 10 -19 J |J. If the rod density pb is too large the rod- rod interaction becomes 
important which might screen the interesting effects |§. Since these interactions scale with 
the second virial coefficient Bi ~ DL 2 a useful dimensionless measure for their strength is 
pi = pbDL 2 . So even if pi must be limited to a small number, 0.1 say, x can in principle be 
made arbitrarily large by choosing large enough aspect ratios L/D. In practice, however, one 
would need L/D ~ O(10 3 ) to obtain x ~ O(10). This is much higher than for the "classical" 



rod-like colloidal particles like the tobacco mosaic virus []T6| but may be achievable with 
microtubules for which D = 25 nm and L can be tens of micrometers. Rather large vesicles 
would have to be used so that the curvature radii are still large compared to L which justifies 
the neglection of higher order terms in the curvature. 

Microtubules and many other mesoscopic rod-like particles are usually polydisperse. Thus 
a useful extension of the present work would be the inclusion of polydispersity which poses 
no fundamental technical problems as long as the interparticle interactions can still be 
neglected. A generalization to soft particle-wall interactions seems to be more difficult as no 
simple analytical expressions for more realistic potentials exist. 

Finally we mention that contrary to what is claimed in Ref. M the corresponding problem 
for disks instead of rods is not completely equivalent. It is easy to convince oneself with a 
coin and a cup that there are configurations where a disk touches the inside of a cylinder 
at two isolated points whose distance is smaller than the disk diameter. The orientational 
constraints due to these configurations obviously cannot be described by replacing the disk 
with an equivalent rod. 
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FIGURES 



FIG. 1. Geometries for the calculation of the available orientational space of a rod near a 
general surface. For positive curvature the rod touches the surface tangentially for small normal 
distances of its center of mass (rod 1) and with its end at larger distances (rod 2). For negative 
curvature only end contact occurs (rod 3) and the rod cannot come closer to the surface than the 
distance z c (rod 4). The value of x max [Eqs. @ and (||)] is the cosine of the angle 6 at contact. 

FIG. 2. Curvature dependence of the surface free energy density w s (ci, C2) for c\ = ccostp and 
C2 = csiwif). The value of c is kept fixed (cL = 0.2) and tp is varied from to 2ir providing a 
continuous path through the four possible sign combinations of c\ and C2 for which four different 
analytical expressions [Eqs. (Pl)- (|l3|) 1 apply. The curve is non-analytic at ip = |, 7r, 2-7T. The 
third derivative with respect to tp diverges on one side of these points. 

FIG. 3. Aspect ratio a and bulging parameter t (see main text) of the equilibrium shape as 
a function of the reduced concentration x of rods inside and outside the vesicle for prolates and 
oblates at two reduced volumina v (co = 0) close to the transition point. 

FIG. 4. Equilibrium shapes for cq = at v = 0.648 which corresponds to the transition point 
if no rods are present. Only one quarter of the contours is drawn, the remaining parts follow by 
symmetry. Note that the symmetry axis is drawn horizontally for the prolates and vertically for 
the oblates. Within each part of the figure the shapes are scaled to the same volume and area. The 
presence of rods outside or inside the vesicle induces a modification of the shape, which, however, 
for the prolates is hardly visible on this scale. 

FIG. 5. Phase boundary between prolates and oblates for cq = as a function of the dimension- 
less rod concentration x outside, inside, or on both sides of the vesicle. Particles outside favor oblate 
shapes, particles inside prolate shapes. If particles are present on both sides their contribution to 
the elastic energy has the Helfrich form, too, so that the equilibrium shapes and the boundary are 
not altered. 

FIG. 6. Phase diagram in the v, cq plane for vesicles in the absence and presence of rods. 
The lines approach the maximum value v = 1 at cq = —6/5 for x = or rods outside, but at 
Co = —6/5(1 — x/64) for rods inside or on both sides (see main text). 

FIG. 7. The same phase diagram as in Fig. || but plotted in terms of the effective spontaneous 
curvature c' (Eq. (|34|)) for rods inside and on both sides of the vesicle. The curve for the latter 
case coincides with that for the free vesicle (x = 0) in this representation. 
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